Observing dynamical SUSY breaking with lattice simulation 
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Abstract. On the basis of the recently developed lattice formulation of supersymmetric theories which keeps a part of the 
supersymmetry, we propose a method of observing dynamical SUSY breaking with lattice simulation. We use Hamiltonian as 
an order parameter and measure the ground state energy as a zero temperature limit of the finite temperature simulation. Our 
method provides a way of obtaining a physical result from the lattice simulation for supersymmetric theories. 
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INTRODUCTION 

In these several years, lattice formulation of supersym- 
metric theories made a great progress IH IH IHH, El 01 
(for the review, see |7]). Not only the formulations but 
the simulations have been akeady started isHoi fToll . 

We have developed a method of observing dynami- 
cal SUSY breaking using a lattice simulation iflll \v^ . 
Since the SUSY is not broken in the perturbation if it is 
not broken in the tree level, a way of observing SUSY 
breaking due to non-perturbative effects is very impor- 
tant. Usually, the Witten index provides such a method. 
But sometimes it is not available. In fact the main tar- 
get in this talk, two-dimensional .y/ = (2,2) super Yang- 
Mills theory (SYM) is such an example. What we know 
is only a conjecture which states it may be broken 1 13]. 

Our method uses the Hamiltonian as an order parame- 
ter. The requirement for the lattice model is an exact Q- 
symmetry on the lattice so it can be widely used in prin- 
ciple. After a very quick explanation of putting SUSY on 
the lattice using Sugino model [|2j], we state our method 
of measuring the Hamiltonian. After confirming that the 
method actually works, we will show the result for the 
super Yang-Mills case. 



LATTICE MODEL 

We use a lattice model with one exactly kept super- 
charge. It seems impossible to put the SUSY on the lat- 
tice, because SUSY algebra contains infinitesimal trans- 
lation but on the lattice we have only finite translations. 
However, what we have realized in the recent develop- 
ment is that it is possible if ^ > 2, especially in low- 
dimensional case. Most of the recent development utilize 
the topological twist. After the twist, we have a scalar su- 
percharge instead of spinors. We can put the scalar on a 
lattice site and keep it exact at finite lattice spacing. Some 



of the simulation have akeady done aiming the check of 
the formulation. 

We adopt such a model proposed by Sugino [2]. The 
target theory in the continuum has twisted ,y\^ — 2 super- 
symmetry. It has the following algebra: 

2 _ ;^ (gauge) q2 _ g^gauge) 



G =5. 



(1) 



{e, go} =2/^0 + 251™, (2) 

where we pick up only a part of the whole algebra, 
g (gauge) (jgjjQjgg infinitesimal gauge transformation with 
the parameter •. After the twist, we have four super- 
charges, a scalar Q, two from a two-dimensional vector 
Qji (M = 0; 1)^ ^nd a pseudo scalar Q. They are nilpotent 
up to gauge transformation. The field components are a 
complex scalar field (0 = 0^), gauge field A^, Majo- 
rana fermions in the twisted basis 77, % ^nd and an 
auxiliary field H. 

The point is that the action has g-exact form S = 
QA, where A is a gauge invariant quantity, and Q can 
be discretized keeping the nilpotency even at a finite 
lattice spacing. Therefore, the lattice action defined as 
g-exact quantity is manifestly g-invariant and enjoys Q- 
symmetry at a finite lattice spacing as long as it is gauge 
invariant. The lattice version of the Q transformation are 
the following: 

QU{x,iJ.) = i\i/^{x)U{x,iJ.), 

QWti{x)=iWyi{x)Wyi 

-i{^{x)-U{x,ii)^{x + fi)U{x,n)-^), 

G0W=o, 

Qxix)^HiX), QHix)^[^{x),x{x)], 
Q^{x)^ri{x), Qri{x)^[^{x),'^{x)], 

where U{x,iJ.) is a gauge link variable. 

Eventually, we have nilpotent Q and g-invariant ac- 
tion. The remaining go, Qi and Q will be automatically 
restored in the continuum limit. 



METHOD 

We use Hamiltonian as the order parameter since it is 
zero if and only if the SUSY is not broken. Therefore, 
it is very sensitive about the choice of the origin of the 
Hamihonian. We use the SUSY algebra to define the 
Hamiltonian. From the algebra Q, we have 



(3) 



where 



r(0) 



is the 0-th component of the Noether current 



for Qo and is a Hamiltonian density. On the lattice we 
have only Q transformation but no Qq transformation. 
Therefore, we discretize the continuum version of the 
Noether current by hand. We know the Q transformation 
so it is straightforward to obtain the Hamiltonian. 

An advantage of measuring the Hamiltonian is that it is 
a one-point function. Compared with two-point functions 
which are used for measuring spectrum, it is numerically 
very cheap and easy to calculate. 

Another point of our method is the boundary condi- 
tion. In the lattice simulation, the lattice size is limited 
and usually we use a periodic condition. For the current 
purpose, however, we need a different choice. 

As usual method for observing spontaneous symmetry 
breaking, we apply an external field conjugate to the 
order parameter. The conjugate to the Hamiltonian is 
the temperature. That is, we impose the anti-periodic 
condition in the time direction for fermion. Therefore we 
break SUSY by boundary condition or equivalently by 
the temperature. Then we take zero temperature limit and 
observe the effect of the breaking is left or not. 

It is interesting to see what would happen if we took 
the periodic conditions. Under the periodic condition, the 
simulation does not work. Since the periodic partition 
function is exactly the Witten index, it is easy to see the 
expectation value of the Hamiltonian is proportional to 
the derivative of the Witten index (exp(— j3//))pBc: 



(H] 



PBC — 



-^(exp(-j3i/))pBc 
(exp(-i3//))pBC 



(4) 



Because the index does not depend on the coupling jSQ, 
the derivative w.r.t. j3 is zero so the expectation value 
of the Hamiltonian seems zero as well. If the SUSY 
is broken, however, the denominator is also zero and 
the expectation value of the Hamiltonian can be non- 
zero value. In the simulation, what we measure is the 
numerator The path integral is replace by an ensemble 
average and the normalization factor is replaced by the 
size of the ensemble. Therefore, the simulation should 
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FIGURE 1. Hamiltonian of SQM versus inverse tempera- 
ture P, in unit a dimensionful parameter m in the potential. The 
plot is after taking the continuum limit. The "not broken" data 
goes to zero as the inverse temperature jS goes large, while the 
"broken" data stays finite. We can easily distinguish these two 
cases so our method actually works. 



give always zero even if the SUSY is broken under the 
periodic boundary condition. 

In fact the expectation value of our g-exact Hamilto- 
nian is always zero under the periodic boundary condi- 
tion, if we ignore the possible zero of the denominator. 
Since the action is g-invariant and the Q is nilpotent, the 
expectation value of Q-exact quantity is zero. This is con- 
sistent with the above argument. 

Eventually, we should measure the ground state en- 
ergy as a zero temperature limit of the Q-exact Hamilto- 
nian. 



RESULT 

As a check of our method, we first investigate Supersym- 
metric Quantum Mechanics (SQM). The known fact is 
that the form of the potential decides whether SUSY is 
broken or not. The lattice model we use was given in |T^] 
and keeps nilpotent Q and Q-exact action. The figure [T] 
shows that our method actually works. Next let us in- 
vestigate the SYM case. See figure|2l We plot results ob- 
tained by two different simulation algorithms, new one 
and old one presented in lilllfl^ . The errors are drasti- 
cally reduced in the new plot and the ground state energy 
seems smallQ Because of the Q-exactness, the periodic 
Hamiltonian should be zero. In fact it is consistent with 
zero within the error (figure O and consistent with the 
discussion based on the Witten index. 



' In the anti-periodic case, fi can be understood as tlie inverse temper- 
ature. 



^ Before giving conclusive statement on the SUSY breaking, we 
should confirm that the lattice model we use actually describes the 
target theory in the continuum H5il - 
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FIGURE 2. Hamiltonian of SYM versus inverse temperature 
jS, in unit of the dimensionful coupling g. The old data was ob- 
tained by reweighting method and the new data was obtained by 
RHMC method. The both indicate small ground state energy. 



CONCLUSION 

We have developed a method of observing dynamical 
SUSY breaking using lattice simulation. Our method 
is available if the lattice model has one exactly kept 
supercharge Q and Q-exact action. We used Hamiltonian 
as the order parameter and measure the ground state 
energy. It actually worked in supersymmetric quantum 
mechanics, a system whether SUSY is broken or not was 
akeady known. Then we applied it to two-dimensional 
super Yang-Mills. It is the first physical application of the 
recent development of the lattice SUSY. The application 
to the other system is straightforward. Now we can use 
the lattice formulation for supersymmetric theories. 
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FIGURE 3. Hamiltonian versus inverse temperature j3, un- 
der the periodic boundary condition. In the reweighting 
method the continuum limit is taken while in the RHMC 
method the lattice spacing is fixed at ag = 0.2357. The values 
should be zero even at the finite lattice spacing. 



Some details of the simulation of SYM are in order. 
The fermion effect is treated by a reweighting method or 
Rational Hybrid Monte Carlo (RHMC). The reweight- 
ing method uses quenched configurations and treats the 
fermion effect as a part of the observable so it is some 
how indirect. The RHMC method uses configurations 
with fermion effects. The simulation codes are developed 
based on the FermiQCD [16] and parameters for rational 
approximation in the RHMC are obtained from We 
fix the physical spacial size to 1 .41 in unit of dimension- 
ful gauge coupling g0The parameters are the following. 
The lattice size is 3 x 6-36 x 12 and the lattice spacing is 
0.0707-0.2357. The number of the independent config- 
urations are 9900-99900 in the reweighting method and 
10-1700 in the RHMC method. 
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^ All the dimensionful quantities are measured in unit of g. 



